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Abstract 

The distribution of magnetic moments in finite ferromagnetic bodies was first in- 
vestigated by Landau and Lifshitz in a famous paper [Phys. Z. Soviet Union, 8, 153 
(1935)], where they obtained the domain structure of a ferromagnetic crystal at low 
temperatures, in the regime of saturated magnetization. In this article, we investigate 
the general properties of ferromagnetic domain walls of uniaxial crystals from the view 
point of the Landau free energy. We present the basic ideas at an introductory level, for 
non-experts. Extending the formalism to the vicinity of the Curie temperature, where 
a general qualitative description by the Landau theory of phase transitions can be ap- 
plied, we find that domain walls tend to suppress the layers, leading to a continuous 
vanishing of the domain structure with anomalous critical exponents. In the saturated 
regime, we discuss the role of domain walls in mesoscopic systems and ferromagnetic 
nanojunctions, relating the observed magnetoresistance with promising applications in 
the recent area of spintronics. 

1 Introduction 

The distribution of the magnetization inside a general ferromagnetic body follows a closed 
flux configuration which leads to the appearance of magnetic domains. For the stripe do- 
main structure, which is common in whiskers, magnetizations of neighboring domains are 
oppositely oriented, separated by 180° Bloch walls of finite width. Inside the domain walls, 
the change of the magnetization is not discrete but smooth. In a pioneering work. Landau 
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and Lifshitz proposed for the first time a tiieory that quantitatively predicted the above con- 
figuration, relating domain sizes and wall width with the dimensions of the body and some 
phenomenological parameters associated with the crystal structure (ill. In their analysis, 
the magnetic energy of the crystal is built as consisting of two terms: one is the exchange 
interaction of the spins, proportional to 

[(Vm,)' + (Vm_v)2 + (Vmj2], (l) 

where m = {mx,mym^) is the magnetic moment, whose absolute value is considered con- 
stant and equal to the saturation value; the other gives the contribution of the magnetic 
anisotropy of the crystal, that competes with the exchange interaction. Assuming an easy 
direction of magnetization along the z-axis, the latter was written as 

P[m2+m2], (2) 

with P > 0. The exchange term Q considered by Landau and Lifshitz can be obtained from 
the classical Heisenberg model, 

N 

(• s 

(with 7 > 0, and 5 indexing the nearest neighbor sites) in the continuous approximation, re- 
placing discrete spin variables S, by a spin density order parameter (magnetization density) 
m(r), spatially averaged in a given local configuration of spins |2|. 

To find the distribution of the magnetization inside the material, one solves a varia- 
tional equation that minimizes the sum of these two contributions for the particular ge- 
ometry and size of the sample, with proper boundary conditions. For the stripe geometry 
and the closed flux configuration, boundary conditions induce the formation of stripe do- 
mains. The solution of the problem leads to a soliton-like pattern for the magnetization 
near domain walls||T|, forming a non-homogeneous phase where the order parameter (mag- 
netization) changes sign alternately from one domain to the next. Due to the ferromagnetic 
interactions, the spins locally tend to be aligned, and one sees that the exchange energy is 
small everywhere, except in the intermediate region between domains, where the magnetic 
moment smoothly changes its orientation to satisfy the anisotropic energy and the bound- 
ary conditions. The overall energy is minimized, since there are no field lines outside the 
sample. The main assumption of this treatment is that the system is ferromagnetically or- 
dered and close to saturation, in the low temperature phase. An extensive description of the 
general properties of magnetic domains and domain walls is available in Ref. 0- 

There are several important applications related to the existence of magnetically ordered 
domain walls, including the recent advances in spintronics for electronic devices. The pur- 
pose of this article is to discuss the structure of the domain walls in a phenomenological 
way, using the simple but efficient tools provided by the Landau analysis on ferromag- 
netism. In the first part of this article we consider the mean field critical regime of domains 
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walls, while the second one is devoted to some effects in the saturated regime, where we 
discuss possible applications. 

To understand how the domain structure behaves in the vicinity of the Curie temper- 
ature, where the magnetic order becomes paramagnetic, we extend the original Landau- 
Lifshitz approach within the more general framework of the Landau theory of phase transi- 
tions. The Landau free energy of a ferromagnetic system is written as a series of the mag- 
netic moment density m = {mx,my,m^), which is considered as the order parameter (OP). 
In the neighborhood of the critical temperature T^, |m| is assumed to be small. This power 
series introduces phenomenological coefficients for the exchange and anisotropy contribu- 
tions, and is written following the general Landau prescription. If one writes the system 
degrees of freedom in terms of the magnetization density m and the magnetic field H, the 
thermodynamics is obtained through the partition function 



exp 



-G(H,r) 

IcbT 



which is written in terms of the Gibbs free energy G(H, T), and is proportional to 

£'(m,H)" 



'Dm) exp 



kBT 



(3) 



(4) 



where ks is the Boltzmann constant and T is the temperature. The symbol Dm means inte- 
gration over all possible configurations of m, with ^(mjH) being the effective Hamiltonian 



£(m,H) = j t/^r/(m,H). 



(5) 



and / the Landau functional. The latter integration is carried for a (i-dimensional system. 
The expression Q yields the partition function as a functional integral of the field. Assum- 
ing that the dimensions of the stripe domains are large in comparison with the lattice con- 
stant (continuous approximation), variations of the magnetization lines inside the material 
are considered smooth and the term that measures the inhomogeneity contribution can be 
treated in first order as a gradient. In other words, we assume that the field m = (m^, my, m^) 
and its derivatives are continuous. In the mean field level, the minimization of the Gibbs 
free energy given by Eq. ^ is done in the saddle point approximation. 



5£fm,Hl 



5m 



0, 



m=(m) 



with the symbol 5 representing the variation of integral (|5Jl with respect to the magnetiza- 
tion, calculated at the saddle point m = (m), where the brackets have the meaning of an 
ensemble average. In this approximation (which neglects the effects of fluctuations around 
the saddle point), the Landau functional /(m,H) plays the role of the free energy density 
and is the relevant quantity to be minimized. In the vicinity of the phase transition, how- 
ever, the fluctuations around the saddle point 8m = m — (m) acquire a major importance for 
physicals systems of low dimension and coordination number Defining ^ as the correlation 
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length, we notice that the mean field analysis breaks down when the amplitude of the fluc- 
tuations measured by the correlation length (5m(r)5m(0)) ^^^'^ is larger than the square 
of the mean field order parameter (m)^ <=< jj4| sj^ce ^ diverges in the critical point, the 
saddle point approximation excludes the region immediately bellow Tc for d <A, where the 
fluctuations predominate. 



1.1 Phenomenological Landau Theory 

It is instructive to develop a few basic ideas regarding the general theory of phase transi- 
tions in a continuous medium. The free energy / is built under symmetry considerations, 
observing the most general form of / which is preserved under all the symmetry operations 
of the OP m that leave the system physically invariant. As we are in the continuum, we 
only need to know the point group symmetry of the OP, which for a ferromagnet depends 
basically on the anisotropy directions in the crystal. In general grounds, the free energy is 
made of three terms. The first one is a series expansion in the OP derivatives K{mi, ^), 
which we call "kinetic" term. 



ij ijk "O ijkl 



where r, are the space coordinates and ni, are the OP components. Not all of these terms are 
to be kept, since first order derivatives on the form /yd'^'r^ are non-extensive when 
integrated in the volume V. In the same way, part of the second term above 



1 



kji) 



drrij dmii- 
hm,-^r— 



^ i^Uk + ^kji) ^ [mm,] 
ijk ^ ^0 



is also non-extensive in the volume. Keeping only the extensive terms, the most general 
form of K up to second order is 13 



1 



ijk 



yrtij 



nik- 



' drj 



ijkl 



ijkl 



drj dri 



(6) 



The first term of Eq. ^ is of special importance in the phenomenological description of 
commensurability transitions |5|. As the free energy of a ferromagnet is invariant by the 
inversion of the coordinate basis in the space, meaning r,- — r,-, the antisymmetric term in 
r,- must be also discarded. The actual form of K is ruled by the physical properties of the 
system and by the point group symmetry of the OP. For a cubic crystal, the tensor Ciju is a 
number and the lowest order term is given by Eq. Q. 

The second term of the free energy / is the potential U (m,), described by a power series 
expansion in terms of the OP components. 



(1) 



L(2) 



+ ■ 



.m; 



0{m 



n,i2, 
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provided that m is small, where ip = 1, runs over the components of the OP. The precise 
form of U also depends on the OP point group. In a ferromagnetic crystal, this term defines 
the energy of crystalline anisotropy. The term of first order is clearly excluded, due to the 
inversion symmetry of the crystal. 

The ferromagnets can be classified according to their axes (or planes) of easy magneti- 
zation. For crystals with one axis of anisotropy, say along the [001] direction (or the z axis), 
the second order term is usually written in the form 



v*^ >{m,.+m^) or — v 



Both forms are equivalent, because they are related by an irrelevant constant v'^^^mfy If 
v^^) > it is said that we have an easy axis of magnetization, while for v^^^ < we have 
an easy plane of magnetization, namely the xy plane. The fourth order invariants are made 
of free combinations of products between the two second order invariants above. We see 
that there are four invariant terms of fourth order entering in the free energy as a linear 
combination. The simplest example of uniaxial ferromagnet corresponds to a crystal with 
tetragonal symmetry. The hexagonal lattice of cobalt is another example of uniaxial ferro- 
magnet, with one easy axis of magnetization perpendicular to the hexagonal lattice. In fact, 
the anisotropy along the hexagonal directions is very small and appears only in the sixth 
order terms [6 |. For biaxial crystals, we need one more free parameter in the second order 
invariant of the xy plane, 

(2) 2 , (2) 2 

In the particular case of isotropic cubic crystals, there are no privileged directions between 
the three principal directions [100], [010] and [001]. In this example, we have three equiv- 
alent easy axes of magnetization. The lowest order term 1/2 = 0. [ml + m^, + mf\ = amg is 
spherically symmetric. In fourth order, the invariants can be written in two different (and 
equivalent) ways 

1 

2' 

related by the constant ^v^^'m^. After these examples, we conclude that we may set the gen- 
eral form of the free energy by simple symmetry arguments only. To be more specific than 
this and specify the values of the parameters left, we need more information from experi- 
ments or from microscopic calculations. At the mean-field level, however, the qualitative 
form of the free energy is in general sufficient for drawing several important conclusions 
about the system. 

The last term of the free energy is due to the inclusion of the external magnetic field 
HQ, 

871 

The expression above satisfies the thermodynamic relation for the total magnetic field B of 
the crystal. 



1/4 = — -v^"^' (m^ + My + m^) or U4 = v^^^ [m^m^ + m^.m^ + niyml) 



^ j/rB,{r) = ^ l^d'rmr) +47i(m,(r))] 



aG(m,H) 



106 



B. Uchoa and G. G. Cabrera 



where the Gibbs potential G has been defined in Eq. (|3ll— (|5ll and M, = Jyd''r{mi) is the 
total magnetization. 

2 Critical Region 

We will concentrate our focus in the simplest case, a crystal with one axis of easy mag- 
netization along the z direction. Anisotropics in the exchange are usually small, and (Q 
is a good approximation, even for axial symmetry. The dominant anisotropy effects in the 
spin Hamiltonian come from the admixture of the spin-orbit coupling into the crystal field. 
Similarly to we will write the second order contribution from the crystalline anisotropy 
in the equivalent form 

[aiml + m^,)+yml], (7) 
with a and y positive constants. In the bulk, far from the domain wall, we have 

rux = My = 0, = ibmo, Vm, = . 

To study the interface criticality, we write the excess free energy relative to a bulk system 
in the following form: 



/(r,r) = a{T) [aiml + mj)+y{ml-ml)] + 

+b [f {ml - m^) + 2ay{ml + m^)ml + o?{ml + m^^f] + (8) 
+c [|Vm;,|2 + |Vmy|2 + |Vm/] . 

As usual, the phenomenological coefficient a(r)for the quadratic term changes sign at the 
critical temperature and is taken as a linear function of T, in the form 

a{T) = (^-^^^ qq = qq t , (9) 

with flo > and t = so is negative in the low temperature phase. We then 

assume a < y, since the z-axis is taken as the easy axis of magnetization. Additionally, we 
will assume that H = 0. 

The above Landau free energy Q is a simple generalization of the commonly used 
expression to study phase transitions (see for instance Ref.|7|, p. 417), with quadratic 
and quartic terms in the order parameter. Due to the magnetic anisotropy, the elementary 
invariants of the axial symmetry group are now (m^ + and m^, and the free energy 
is built from them|5i|. The fourth order term in the free energy can be interpreted as the 
contribution of the spin quadrupole interaction. To minimize the number of parameters, 
the quartic term is written as the square of the quadratic one. So we are left with the set 
{a{T),b,c,a,y) of free parameters, where a and y are determined by the crystal field of 
the ferromagnet. The others are usual Landau parameters, with b and c slowly varying 
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with temperature, even at the critical point. We are assuming that they are constant. The 
behavior at large |m| is dominated by the quartic term in (ISj. This requires the constant b 
to be positive, in order to have minima of the free energy at m 7^0, for the low temperature 
phase (where a < 0). The constant c is also positive, since it costs some energy to create an 
interface. 

The original treatment of Landau and Lif shitz argues that the minimization of the energy 
has to be done in two different regions of the crystal: i) in the intermediate region between 
domains (interface), where the contribution of the inhomogeneity cannot be neglected; and 
ii) in the region close to the surfaces, where the closed flux configuration of the field requires 
appropriate boundary conditions. We will use here the same argument, but for the more 
general free energy given by (HJ. 

2.1 Domain Wall Magnetization 

We shall concentrate first on the interface region. Without loss of generality, we assume that 
the crystal is infinite and the interface between domains is in the 3'z-plane. In the absence 
of a magnetic field, due to symmetry, the spins are all in the 3'z-plane. Far from the wall, 
they are aligned with the z-axis (parallel or anti-parallel). Close to the walls, we assume 
that they are deflected by an angle 9 from the z-axis, with components: 

mx = 0, my = mosin6, m^=mQCO&Q, (10) 

where 6 = 6(jc) is a function of x only. Note that wiq = |m| is temperature dependent, but 
considered uniform in space. The substitution of (flOl into ^ leads to 

/ [Q{x),Q'{xy,T]=cml{Q'f + a{T)[a+{y-a)cos^d]ml-a{T)yml+ 

(11) 

+b ['f cos"^ 6 + 2Yasin^ 6cos^ 6 + sin"^ 6]/Mq - tfrn^. 

In order to apply the Euler variational principle to the action one notes that the 'La- 
grangian ' / [6(x),6'(x);r] does not explicitly depend on the variable x. In this case, the 
Euler-Lagrange equation can be written as 

d 



dx 

which reduces to 







cmo(e')^ = K + a{T)[a+{y-a)cos^Q]ml-a{T)yml + 

(12) 

+b [-f cos"^ e + lyasin^ Ocos^ 6 + sin"* Q]mQ - b-frriQ. 

where K is a. constant to be evaluated using the boundary conditions (BC). Far from the 
interface between domains, we set the spins asymptotically aligned with the z-axis, remem- 
bering that z is our easy magnetization direction. To satisfy the closure configuration, we 
impose 
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(13) 



and 



e' = , ^ ±00 , 
which result inK = 0. The equation to be solved now is 



(14) 



cmliQ'f =a{T)[{a-y)sm^Q]ml + bm'^ [2Yasin2ecos^e + a^sin'*e + y^(cos'^e- 1)] , 

(15) 

which means that the free energy is minimized when the exchange energy density is equal 
to the anisotropy one. Equation (031 can be written in the form 



(6')^ = sin2e (A + Bcos^e) , 



whose solution can be given as 



where A and B are 



1 + f cose 
1 + f cos2e 



-tanh[WA + S] , 



(16) 



(17) 



B 



(a-Y) 



1 H m^{a + i) 



a-Y) 



-n?{a-i) 



(18) 



The above equations give the distribution of the magnetic moment density in a general 
crystal between two neighboring domains. There are two equivalent solutions with opposite 
helicities. The width of the interface (Bloch wall) is given by 



1 



Va+b' 



(19) 



To check the consistency of (fTTl . we remind the reader that a{T) is a negative scalar 
function for T < Tc. Noting that a < y (easy z— direction), the condition (A + B) > 0, will 
impose limitations for the lower bound of a. In our discussion, we are considering the 
limit of small anisotropy, for which a < Y we will show below, mean field implies 
that bm^/a{T) is finite at the critical point, which determines the critical exponent of the 
magnetization. We will discuss those points later on. 
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The symmetry of the magnetization planes in the presence of an uniaxial anisotropy is 
represented by the C21, point group, which admits two invariant terms, and /MqCos26 in 
the expansion of the free energy |5 1. This is easily noticed by decomposing the anisotropy 
term U2, (asin^ + ycos^ 9) [see (0 and (flOb l. into the polynomial basis formed by the 
two invariants, giving 

^ [(a + y)mg + (Y-a)mgcos(2e)] . 

Clearly, the full symmetry term htq due to the paramagnetic phase (0(2) point group) com- 
petes with the other term representing the symmetry of the domain ordered phase (Cjv point 
group). The stability of the domain walls is ruled by the anisotropy parameter £, whose 
proper expression is 

y- a 



0<8: 



y+a 



< 1. 



Assuming that the anisotropy is small (e <^ 1), the ratio 

(A + B) >0, 

B labnin e 
A 



is also small provided that 



and the solution can be expanded in terms of j 



1 + 



2A, 



j COS 6 


B . 


1 — cos^e 




2A 



■tanh[v^(l + — );c] 



(20) 



Fig. 1 shows the distribution of the magnetization angle 6 through the interface. The 
quantity (A + B) vanishes at the critical point yielding the limit X —>■ 00 and the solutions 
9 = ±1 everywhere. At first sight, this solution may seem inconsistent with the boundary 
condition imposed. However, we observe that the magnetization intensity goes to zero at T^. 
In any case, this behavior signals that something odd is happening at the critical point and 
questions the validity of mean field solutions there. The limit X — > 00 means that magnetic 
fluctuations are paramount at Tc. 



2.2 Surface Energy 

Next we consider a finite crystal, where the domain structure is organized in layers, as 
shown in Fig.2. We already know the spin distribution through the Bloch wall, and want 
to calculate the width of magnetic domains in a finite volume. This is done in a variational 
way, as in the original contribution by Landau and LifshitzQ. We note that the flux closure 
condition induces the formation of small domains near the surfaces, where the magnetiza- 
tion points perpendicular to the easy direction. 

We proceed to the calculation of the wall energy. If I, and ly are the dimensions of 
the crystal in the z, x and y directions respectively, the energy associated with one interface 
(wall) between two domains is 
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-2-10123 
X (length unit) 



Figure 1: Magnetization profile in the vicinity of a domain wall. The solid curve represents 
the solution in the saturated regime at low temperature (asymptotic case for B ^ and 
A = 8, now normalized to adimensional parameters). The dashed and dot-dashed curves 
were obtained solving d20b . keeping the anisotropy fixed and varying the temperature in 
direction to Tc, where both, A and B vanish (dashed for B = 0.1, A = land dot-dashed for 
B = 0.05, A = 0.5, respectively). 



/oo 
<ix/(m,T) = 
-oo 

/oo 
dx [cml{Q'f + a{T) [a+{y- a) cos^ e]mg - a{T)yml+ 
-oo 

+b [y^ cos"* e + lya^in^ Ocos^ 6 sin"^ 6] nig - bfmQ\ (21) 

where the limits of integration have been extended to ]— oo,oo[ , considering that the layer 
width is much larger that the wall region. The dominant contribution to this integral is 
concentrated inside domain walls, where 6 is close to 7i/2. This way, we may neglect the 
cubic term in cos(6) in equation d20b . Substituting expressions (fTSl into ^T\^ yields an 
expression only in terms of the anisotropy 

/oo 
dx sin^O \a{T)m^+ 
-oo 

-hZ7[(a-FY) - (a-Y)cos^e]mo] , 

where we have used the form given by (fT6b . To lowest order in the anisotropy 8 we get the 
result 

/°° dx 
7 : 
~cosh^ [x^/A 
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It then follows that the total energy of {Ix/d) walls is 



-c m. 



where d is the layer width. 

Now, we turn to the surfaces for the flux closure configuration. The ferromagnetic phase 
of the body produces some amount of field given by the Maxwell equation V • (H + 47im) = 
inside the sample, and V • H = outside. Near the surface and far from domain walls, 
the exchange contribution to the energy is zero. The closed flux configuration, with no 
field lines outside the sample, is favored when the anisotropy is small. Under the above 
conditions, we get 

H = 0, and V • m = , 

what means that the spin distribution satisfies the boundary condition m • n = 0, with n 
normal to the surface. This way, surface poles are avoided and the global spin field inside the 
crystal has no singular points as shown in Fig. 2. For surface domains, the main contribution 
to the energy density comes from the anisotropy. At the surface, we get 



Mr = iiin 



0, 



niy =0 



and going back to Q we get the energy density 

/surface = a{T)ml{a - j) 



1 + — ^(a + Y) 



which is again proportional to the anisotropy £. For a finite sample, we have two opposite 
surfaces at z = and z = —I- Associating the volume {lyd^ /4) to a single surface domain, 
and summing over {2lx/d) of such domains, the total surface energy is given by 

1 2 

El = -Uyd MqCA . 

We then minimize the total energy E = E\+E2 simultaneously in relation to d and mo 



dd 



0, 



=0 

3mo 



yielding the results 



d = 2V2lJJ\, 



bniQ 

I 4f) 



(22) 



4 1 



7 (a+Y) 
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Figure 2: Magnetic permeability distribution in the crystal for flux closure. 



The quantity A, defined in (fT^ . vanishes at the critical point, thus causing the divergence 
of the domain width at Tc. The associated critical exponent is different from the one that 
gives the divergence of the wall width X in ( fWl 



In the expression above we encounter the quantity 




D = \ + 



2ybml 



that has to be positive in order to get real solutions. From (1^ . we see that this is the case 
for y/7 < a < y. So, we are on safe grounds for small anisotropy, a < y. We note that the 
critical exponent for the magnetization is the same as in the homogeneous case (P = 1 /2), 
but the full expression is different from the one obtained in conventional mean field theory. 
It is worth to note that the isotropic, homogeneous case, satisfies D = 0, i.e. no domains are 
present. 

From d22b . we get the result 







mo 



7fc(a+Y) 



,T>Tc 

,T<Ta 



The two scales, d and X, diverge at the critical point, but the divergence of the wall width is 
faster (exponent —1/2) than the one for the domain width (exponent —1/4), thus showing 
that the critical point marks the onset of strong magnetic fluctuations. 

For magnetic thin films, the Bloch walls produce magnetostatic poles in the surface 
of the film, enhancing enormously the energy of the crystal when the thickness is smaller 
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than ~ lO^^cm [8|. In general, the crystal prefers to deflect the spins of the wall along 
the directions parallel to the thin film surface, giving rise to a Neel wall, when the easy 
magnetization direction is in the plane of the film, as shown schematically in Fig. 3. 

Despite it is tempting to generalize our previous conclusions to this system, the long- 
range dipolar interactions (which we have neglected so far) are relevant in an infinite slab 
of small thickness, not mattering how weak the dipolar interaction is in comparison to 
the nearest neighbor exchange interaction f^. To illustrate how important this interaction 
might be, in some ferromagnets, the dipolar interactions are responsible for the appearance 
of domains in the complete absence of crystalline field effects, when the free energy is fully 
isotropic in a given plane of magnetization, for example. In these systems, the exchange 
competes with the dipolar interactions, stabilizing a non-homogeneous phase with a finite 
wave-vector in an arbitrary direction of the magnetization plane. Because of the large phase 
space of all the degenerate wave- vector directions allowed by symmetry, the fluctuations are 
strongly enhanced at Tc, giving rise to the fluctuation induced first order transitions studied 
by Brazovskii 1101 . In this specific case, the study of the phase transition is totally out 
of the scope of the mean field analysis. Here, we have focused our study in ordinary 3D 
ferromagnets, where the dipolar contribution is not a relevant interaction. 



£ 
c 

-1 



^ < w ^ 

< ^ *K 

^ < w ^ 

^ < w ^ 

^ < w ^ 





x/A 



Figure 3: Illustration of a 180° Neel wall. 



2.3 Critical Exponents 

When temperature approaches Tc from below, a{T) as defined in goes to zero linearly, 
and the entire system is affected by long range fluctuations. At Tc, the correlation length 
diverges and the system is scale invariant|7 1. In our model, this is signaled by the divergence 
of the wall width X, which plays the role of the length scale and goes to infinity under 
the power law X°<^ kl -■ Since the domain width behaves as d <x \t\^^, we find that the 
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wall region enlarges and 'compresses' the domains as temperature raises in direction to Tc, 
meaning that domains firstly lose their identity and finally disappear at the critical point. 
We observe that all the critical exponents in our treatment refer to behaviors below Tc (in 
critical phenomena, one distinguishes the behaviors above and below ). 

In what follows, we advance conclusions derived from the previous calculation. To 
calculate the heat capacity, we compute the total internal energy of the system as 

E = 2cV2llJymlA'i, 

once the minimization process is done. Note that this is the excess energy relative to 
the bulk system, associated with surfaces and interfaces. From this, we calculate the spe- 
cific heat C = and find that C oc \t\^ diverges at f = 0, with the critical exponent 

a = 1/4, which is different from the standard mean field approximation = for homo- 
geneous systems (bulk magnetization) and also different from standard mean field results 
for interfaces = 1 /2)yi |. The different critical behavior comes here from the criticality 
of the domain width d. 

We also know from experiments that the static magnetic susceptibility diverges at the 
critical point. In our calculation we get 

1=[—^] oc |?| 3, for?<l , 

and therefore Y= |, which also disagrees with the standard mean field value y= 1. Applying 
the Josephson, Fisher and Widom scaling lawsf71 fT2ll 

vd = 2-a and y = v(2-ri) = p(5- 1) 

for the correlation length ^ oc t^^ , for the equation of state M = Jydf^x{m) oc H^^, and for 
the correlation function (F) power law p = d — 2 + r\, with 

we find that V = rj = f and 8 = |. 

The different exponents and the singularity in the heat capacity are explained by the 
finite size effects manifested in the geometry of the ferromagnet surfaces, which lead to the 
formation of Bloch walls between magnetic domains. This statement can be easily tested 
by removing the closed flux boundary conditions in the thermodynamic limit. In this limit, 
the magnetization will be homogeneously distributed as if the crystal had a single domain 
with the size of the system, i.e. d ~ 4 not mattering the temperature. We regain all the 
standard mean field exponents in this situation. 

3 Mesoscopic Systems 

The main effect resulting from the introduction of an external magnetic field H in the fer- 
romagnet is the displacement of the domain walls from the zero field equiUbrium position. 
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As the bulk energy of each domain is enhanced or reduced by the magnetic coupling of the 
magnetization with the magnetic field, 



depending on the orientation of the spins with respect to H, the domain walls are dislocated 
in the direction that reduces the size of energetically unfavorable domains, enlarging the 
favorable ones and reestablishing the equilibrium. As shown by Landau and Lifshitz 1 1 1, the 
dynamical interaction of the magnetization with the magnetic field is driven by the Zeeman 
effect, where the spins start to precess around H as free moments, and by the relativistic 
interaction (with coupling rj), following the equation of motion 



where ^a = e/ mc, with e the electron charge, m the electronic mass and c the speed of light. 
The dot represents a time derivative. Considering the case where H is oriented along the 
easy magnetization direction, they found the domain walls speed of displacement to be 



employing the notation of sec. 2 for the saturated regime (a = const, and b = 0). If the 
magnetic field is constant in time, the velocity will be reduced to zero by the magnetic 
pressure of the bulk in a second moment, as the system approaches the equilibrium. On the 
other hand, if we orient the external magnetic field transversely to the easy magnetization 
axis, the domain walls will remain in equilibrium with the magnetization of the bulk. The 
application of external magnetic fields on domain walls may have important applications in 
the fabrication of switches for bulk spin polarized electrons in mesoscopic systems. 

The scattering of an electron through a 180° domain wall barrier comprehends two 
cases of technological interest for the fabrication of switches. Consider for example one 
electron of spin up, aligned with the magnetization of the bulk, flowing in direction to a 
domain wall. If the time of flight of the electron through the wall is long enough, its spin 
will be adiabatically deflected and the scattered electron will be transmitted to the other side 
of the barrier with the spin flipped down, causing no additional cost of energy to the bulk. 
On the other hand, if the time of flight is too short, (considering that the electron is too fast 
or the barrier is too narrow) the spin will not have time enough to be deflected by the spiral 
of magnetization and the electron (in the classical picture) will be reflected back from the 
barrier, in order to save energy from the bulk. 

As shown by Cabrera and Falicov fT3], the domain walls may respond for a large change 
in the resistivity of metallic ferromagnets when the difference in the density of states be- 
tween majority and minority spins at the Fermi surface is also large. Consider that we have 
two spherical Fermi surfaces of different sizes, for example, one for spin up electrons and 



-m H 





116 



B. Uchoa and G. G. Cabrera 



another for spin down ones. In the situation where the radius of one of the bulk Fermi sur- 
faces is very small (minority spin) and the other very large (majority one), the majority spin 
electrons in one side of the barrier will find few channels to tunnel ballistically to the other 
side, where they would occupy the minority spin states. This way, if the wall is sufficiently 
narrow in comparison to the electron mean free path, the electrons will most probably be 
backscattered by the barrier. The change in the bulk resistivity predicted, however, is not 
so large as to induce giant and colossal magnetoresistances. In contrast, the magnetoresis- 
tance can grow several orders of magnitude when the width of the walls becomes small in 
comparison to its size in bulk, as it has been observed experimentally in magnetic nano- 
j unctions. 

A very interesting physics shows up in mesoscopic systems when a domain wall is 
geometrically constrained by a small constriction separating two ferromagnetic bulks of 
wider cross section. Bruno |14| derived the somehow remarkable result that when the 
constriction cross section is much smaller than the bulk one, the domain walls are practically 
independent on the specific characteristics of the material, including crystalline anisotropy 
and exchange stiffness, and depend only on the geometry of the constriction. The free 
energy of the system is in the form 



f[Q{x),Q'{x)] = [cml{Q'f + U{Q)]S{x), 



(23) 



where S{x) defines the geometry of the constriction and U is the crystalline anisotropy term. 
Proceeding with the minimization, one finds through the usual Euler-Lagrange equation 



d_dl 
dxdB' 



ae 







that 



e" + e 



1 du 

2cml 36 



0. 



Bruno proposes that if the second term is much larger than the third and in addition if S is 
integrable, meaning if J^^S{x)dx is finite, then we may drop the crystalline contribution. 
Defining the wall width according to the criterion 
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he finds the simple results: 



e(;c) = 71 



J-^S-^{x)dx~2 



and 
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which are defined by the constriction only. In the most general case, where the crystalline 
term U is present, he finds for different models of constriction that the domain wall is still 
driven by S, and not by U or the exchange stiffness c. 

When a magnetic field is applied along the easy magnetization axis, the domain wall 
inside the constriction suffers a finite displacement proportional to the applied field, where 
the new position of equilibrium is achieved thanks to the magnetostatic pressure played 
by the constriction geometry fl?ll . The localization of domain walls in nanojuctions are a 
theoretical possibility for the fabrication of fast switches demanded for spintronic devices. 

In a wide sense, the word spintronics means manipulation and control of the spin de- 
grees of freedom in condensed matter systems, and to use this knowledge to make useful 
devices fT6l . Here, we focus on spin transport properties in ferromagnetic metals. We have 
already discussed the magnetoresistance effect (MR) due to the domain structure and the 
scattering of electrons from Bloch walls. The resistance of a ferromagnetic sample changes 
with the application of a small magnetic field. The field reorients the domains, and for a 
large enough value, it saturates the sample. The saturated configuration has a smaller resis- 
tance (negative magnetoresistance) than the one with domains, where the magnetization is 
not homogeneous (for a 180" wall, the magnetization of one domain is antiparallel to the 
magnetization of the adjacent domain)|13|. This simple system suggests the invention of 
"artificial" devices where one could monitor the magnetic configuration through the appli- 
cation of small magnetic fields. Based on the same physics, magnetic tunneling junctions 
(MTJ)|17| and metallic multilayers (MML) structures 1 18 1 were conceived, in order to op- 
erate as spin valves, due to the 'giant' magnetoresistance (GMR)effect displayed by those 
systems. ^ In a typical MTJ, two ferromagnetic metallic electrodes are separated by a narrow 
nonmagnetic insulating layer For MML, two or more ferromagnetic layers are separated 
by nonmagnetic metallic spacers. In both structures, one can pin the magnetization of one 
of the electrodes using an additional magnetic layer which is strongly exchanged-biased 
with the ferromagnet llQI . The other ferromagnetic electrode is left 'free' to orient its mag- 
netization with small applied fields. One can then change the resistance of the device by 
manipulating the relative orientation of the magnetization of both electrodes, phenomenon 
that has been named as spin-valve effect and is currently used in commercial applications 
for magnetic recording 120 1 l2Tl . For technological applications, ideal devices should have 
a GMR as large as possible, obtained at room temperature with small applied fields. Ex- 
perimental measurements of GMR in magnetic multilayers are usually performed for the 
current-in-plane (CIP) geometry, where the electric current flows parallel to the layers. For 
this setup, the resistance is fairly large and can be measured using standard techniques. The 
other possibility, much more complicated for experimental implementation, is to do the 
measurement with the current perpendicular to the layer plane, the so called CPP geometry 
(CPP stands for current-perpendicular-to-the-plane). In this case, the resistance is very low 

^The MR effect is called GMR for values exceding 10%, which are one order of magnitude larger than the 
typical values of the anisotropic magnetoresistance used commercially. MTJ prototypes currently produced 
nowadays present a GMR in the range 25%-30% at room temperature. In MML systems, the MR may be as 
large as 50% - 65% , but saturating fields are substantially higher (about 1 T in the original system of Ref. 
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and can only be measured using extremely sensitive techniques 1 22 1. 

Domain walls are thought to play a relevant role in recent GMR experiments in mag- 
netic nanocontacts. Large values of MR, of the order of 300%-3000% at room temperature 
and for fields of about 100 Oe, have been reported in the literature I23l 1241 . The MR val- 
ues found are attributed to strong electron scattering from narrow domain walls which are 
formed in the contact region. Due to the constricted geometry and the rapid variation of 
the magnetization across the domain wall, the electron spin cannot follow adiabatically the 
local magnetization, as it is the case in bulk f erromagnets II 1 3 1 . 

4 Conclusion 

We have developed an illustrative application of the Landau mean field theory for phase 
transitions in ferromagnetic bodies. The advantage of this procedure resides in the phe- 
nomenological nature of the free energy parameters, allowing us to calculate equations of 
state without stating a specific microscopic model. The classical Landau approach can then 
be extended to describe the formation of domains, with the corresponding change of the 
critical behavior, still within the mean field approximation. Despite this theory is in fact 
very crude to describe the neighborhood of the critical point for systems whose specific 
heat diverges at the phase transition [12|, the theory is also known for making remarkable 
qualitative predictions, with deep insights on the physics of critical phenomena, even in 
situations that extrapolate its region of validity|25 |. In particular, the theory provides its 
own criterion of failure and suggests new procedures to correctly describe the physics at the 
critical point, where fluctuations must be included|26 1. 

We have shown that the problem of domains in ferromagnets is driven by the presence of 
two different length scales competing in the vicinity of the critical point. On one side, there 
is the bulk correlation length which drives the system to a classical mean field ferromagnetic 
behavior, and on the other, the length scale of the interfaces between domains. The competi- 
tion between them produces two major effects: i) a new critical behavior, where the specific 
heat is strongly enhanced near Tc from below in comparison to the bulk homogeneous limit, 
but is still small in comparison to the limit of surfaces, giving an intermediary critical expo- 
nent ttfi < a < a^; and ii) a raise in the magnetoresistance in the ordered phase, specially 
for ferromagnets with a layered domain wall geometry. Despite the study of magnetic do- 
mains is a rather old subject, it has recently opened a very promising field in technology. 
The transport properties of ferromagnetic domains suggest several interesting appUcations 
in spintronic devices, which we have briefly discussed. 
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